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Abstract 



> 

' We introduce several coordinate invariant statistical procedures in 

ly-j I order to test for local alignment of polarizations. A large scale align- 

ment of optical polarizations from distant QSOs has recently been 
I observed by Hutsemekers and collaborators. The new statistical pro- 

' cedures are based on comparing polarizations at different angular coor- 

dinates by making a parallel transport. The results of these statistical 
• procedures continue to support the existence of the large scale align- 

ment effect in the QSO optical polarization data. The alignment is 
Q ' found to be much more pronounced in the data sample with low de- 

^ I grees of polarization p < 2%. This suggests that the alignment may be 

^ • attributed to some propagation effect. The distance scale over which 

' the alignment effect is dominant is found to be of order 1 Gpc. We also 

find that a very large scale alignment is present in the large redshift, 
z > 1, data sample. Infact the data sample with z > 1 appears to be 
}J] ■ aligned over the entire celestial sphere. We discuss possible physical 

I effects, such as extinction and pseudoscalar-photon mixing, which may 

be responsible for the observations. 

Key words : polarization: magnetic fields: elementary par- 
ticles: methods - data analysis, statistical: quasars - general. 



1 Introduction 

The optical polarizations from QSOs appear to be aligned with one an- 
other over very large distances (Hutsemekers 1998). This effect was further 
confirmed by Hutsemekers &: Lamy (2001) using a larger data set. It was 
observed that the polarizations from QSOs in any particular spatial region 
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or a patch have a tendency to ahgn with one another, without any evidence 
of a large scale anisotropy. The effect was found to be rcdshift dependent, 
namely the patches which are aligned are delimited both in angular as well 
as radial (redshift) coordinates. A very striking alignment was found in the 
region, called Al by Hutsemekers (1998), delimited in Right Ascension by 
11^15™ < RA < 14'^29™ and in redshift by 1.0 < z < 2.3. By using several 
statistical tests (Hutsemekers 1998; Hutsemekers &: Lamy 2001) the authors 
were able to rule out the hypothesis of uniform distribution at approximately 
0.1% significance level. The basic idea behind these statistical procedures 
is to test for the dispersion in the orientations of polarizations in any small 
neighbourhood. The tests require comparing the polarization vectors lo- 
cated at different angular coordinates. These polarizations are specified by 
the angles they make with respect to the local meridians (Hutsemekers 1998). 
The statistical tests used by Hutsemekers (1998) are, however, dependent 
on the precise position of the pole used for defining the coordinate system 
and hence it is difficult to interpret the results obtained by these tests. In 
the present paper we introduce some statistical tests which are independent 
of the coordinate system and apply them to the QSO optical polarization 
data. 

2 Parallel Transport 

The basic problem can be explained as follows. Consider two tangent vectors 
Vi and 1)2 located at two different points Pi and P2 on the surface of the 
sphere. These vectors are specified at each point on the sphere in terms of 
the angle they make with respect to the local unit vector (f>. Here we have 
chosen the spherical polar coordinate system defined by the unit vectors 
{f,9,(/)). Let ai be the angle made by the vector with respect to the 
local unit vector (pi at the location Pi. Similarly let 02 be the angle made 
by the vector V2 with respect to the local unit vector ^2 at the location P2. 
The two vectors in their local coordinates are given by 

Vi = cos aicpi + sin ai^i 
V2 = cos a2(j)2 + sin 02^2 

In order to compare these two vectors we may use the usual dot product, 

Vi • V2 = cos(ai) cos(a2) + sin(ai) sin(Q;2) . (1) 

However this dot product is not invariant under coordinate transformation 
and hence a statistical procedure based on this will be coordinate dependent. 
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A coordinate invariant dot product may be obtained if we parallel transport 
V\ to P2 and then take the dot product of the transported vector vi' with 
V2. Let the parallel transported vector v^' make the angle a'l with respect 
to 4>2- We are interested in finding a'l which can then allow us to calculate 
the angle a'^ — 0:2 between the parallel transported vector v-l' and the vector 

Let fi and r2 be unit radial vectors at the points Pi and P2 respectively. 
The unit vector s perpendicular to the plane containing these two radial 
vectors is given by 

5 = A^ (2) 

Therefore the unit tangent vectors ti and ^2 at the two points Pi and P2 
respectively along the great circle passing through these two points are given 

ii = s X ri (3) 
^2 = s X r2 (4) 
In terms of the local basis {6, (f) these two vectors are given by 

ii = Bi • iJi + (^1 • fi0i (5) 
(2 = 02 ■ ^2 + h ■ (6) 

where 

- -V — sin^i cos ^2 + cos 01 sin ^2 cos(^2 — 
01-^ = /-, ^ 

VI - (n • r2Y 

~ sin6l2sin((/>2 - 0i) , , 

a/I - (n • r2Y 

- — sin02 cos^i + COS02 sin^i cos((?!)2 — (/)i) 

02-t2 = ===== (9) 

VI - [ri ■ r2Y 

<p2-t2 = j===== (iUj 

VI - (n • r2Y 

As we parallel transport a vector along the great circle its angle with 
respect to the tangent to the great circle remains fixed. Hence in order to 
determine the angle by which the vector has turned due to parallel transport 
we only need to find the orientation of £1 and £2 with respect to the local 
basis at the points Pi and P2 respectively. The angle ^1 between i\ and 
<l)\ is given by cos~^{(f)i ■ ti). This gives ^1 upto an overall addition of tt, 
which is fixed by the sign of 6\ ■ ti. If 6*1 • ti < then ^1 lies in the third or 
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Figure 1: The curve is the great circle passing through points Pi and P2. ti 
and t2 are the tangent vectors at points Pi and P2 respectively. Polarization 
vector vi from location Pi is parallel transported to location P2. Parallel 
transported vector is indicated by v[. 
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fourth quadrand. The angle ^2 is obtained in a similar manner. Hence the 
transported vector v-^' makes an angle (cti +^2 — ^i) with respect to the 4>2- 
Therefore the rule for comparing vectors at two different locations on the 
sphere is given by the generalized dot product, 

ViQv2 = \vi\\v2\ cos(q;i — a2 + ^2 — ^1) , (H) 

i.e. the dotproduct between the parallel transported vector vi and the 
vector V2- 

In our discussion so far we have explained the procedure for comparing 
the vectors at two different points. The polarization orientations ai and 02 
at two different points Pi and P2 respectively can also be compared in a 
similar manner. We associate the unit vectors t^i = [cos(2ai), sin(2ai)] and 
V2 = [cos(2a2), sin(2a2)] with the polarizations at Pi and P2 respectively. 
We parallel transport the polarization cti to the location P2. The angle 
between the parallel transported polarization and the unit vector ^2 along 
the local (p axis at P2 is equal to a'l = (ai + — We can assign a unit 
vector i)^ = [cos(2a'^), sin(2a'^)] with this transported polarization. Wc then 
define the generalized product (ai, 02) between the two polarizations as the 
ordinary dot product between the two vector and V2 

{ai,a2) = v[ v2 = cos[2(ai - 02 + ^ - Ci)] • (12) 



3 Statistical Tests for Alignment 

The simplest statistics to test the alignment of vectors is to calculate the 
dispersion of the polarization with respect to its nearest neighbours. Con- 
sider the riy nearest neighbours of the vector Vi = [cos(20i), sin(2^j)] located 
at the i-th location, including the vector Vi itself. Here 9i is the polariza- 
tion orientation. The orientation angle may be measured with respect to 
the local meridian or the local latitude. Our final results depend only on 
the differences of angles and hence our formulae are directly applicable for 
either definition. The nearest neighbours are determined by computing the 
relative distance in three dimensions with the co- moving radial distance r{z) 
given in terms of the redshift z by the standard relation, 

r(z) = ^ fl-(l + z)-V2l (13) 
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where Hq is the Hubble constant and c is the velocity of hght. A measure 
of the dispersion di of these set of vectors is given by, 



k=l 



5] cos [20 - 2 (0,, + Afc^,)] 



(14) 



fc=i 



where A^^j is the angle by which the polarization orientation angle 
changes after being parallel transported to the position of the polarization 
angle 9i and is a measure of the mean polarization at the i-th position. 
The symbol {0,6k) is defined in Eq. 1121 Hence di is calculated by parallel 
transporting the nearest neighbours of the polarization 9i to the i-th 
position and then taking the dot product of the resulting vectors v'j^ = 
[cos(20fc + 2Afc_+j)] with the vector v{6) = [cos(20), sin(20)]. The nearest 
neighbours include the vector Vi itself. The magnitude of di{9) is then 
maximized with respect to 9. The resulting vector v{9) gives a measure of 
the local mean direction and the magnitude of the maximized di{0) gives 
the measure of dispersion. We define the statistic as follows, 



where the sum is over the entire data sample. A large value of di implies 
small dispersion and hence a large Sjy would imply a strong alignment be- 
tween the polarization vectors. 

It is clear that the mean value of this statistic for a random sample is 
proportional to 1/y^n^. We verify this explicitly by numerical simulations. 
For riv = 25 we find that the mean value = 0.178. In Fig. [21 we show 
the histogram of the statistic obtained by using 10000 random samples 
with = 25. The distribution tends to a normal distribution for large 
data samples. The distribution for different values of Uy was found to be 
identical except for a shift in the mean position, which as mentioned above, 
varies as l/-yn^. In our calculations we evaluated the statistical significance 
level (S. L.) for each case by explicit numerical simulation using a large 
number of random samples. The S. L. is defined as the probability that the 
value of the statistic obtained from the data may be obtained as a statistical 
fluctuation from a random sample. It is numerically evaluated by computing 
the statistic (S*) of the data sample with the statistic S obtained from a 
large number of random samples. The S. L. is a probability that S < S*. 




i=l 



(15) 
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Figure 2: The histogram of the statistic 5^ using 10000 random sample 
with the number of nearest neighbours = 25. The mean position of the 
statistic for a random sample is shown by the verticle line. 
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Our second statistical test is the modified Andrews & Wasserman test, 
which was used by Hutsemekers (1998) and Hutsemekers & Lamy (2001). 
We first review this test. For each vector Vj = (cos 20j , sin 20j ) the mean 
resultant vector with respect to its nearest neighbours can be written as 

Vj = ^{f^cos2ek,Y,sm2eA . (16) 

^fc=l k=l I 



where is a normalization constant such that the vector Vj is normalized to 
unity. The sum over the nearest neighbours Uy excludes the j-th vector itself. 
We then define the mean direction 9j such that Vj = (cos 29 j, sin 26 j). A 
measure of the alignment of vector Vi = (cos 20j, sin 29i) with the n„ nearest 
neighbours of the vector vj = (cos 29 j, sin 29 j) is given by the dot product 
Dij = Vi ■ Vj. In order to evaluate the statistic, the Dij, j = 1, ..,n, values 
are sorted in the ascending order for each i. The rank TZi of Dij=i is then 
evaluated. If n is the number of sources in the sample then the statistic is 
given as 

^ ^J_^'R^-jn + iy2 

As shown by Bietenholz (1986), the statistic Zc approximately follows a nor- 
mal distribution. However we did not use this asymptotic result in obtaining 
our S. L. but evaluated it explicitly by using a large number of random sam- 
ples. The Andrews and Wasserman test can be modified (Hutsemekers 1998) 
by using the unnormalized vector Vj, i.e. by setting the normalization con- 
stant M to unity. The resulting statistic is denoted by (Hutsemekers 
1998). 

As explained earlier this statistic is not coordinate invariant since it com- 
pares the vectors at different locations without making a parallel transport. 
We can modify this statistic by using the generalized product, Eq. ^1 to 
compare polarizations at different positions. We therefore express Dij di- 
rectly in terms of the coordinate invariant quantity {9i,9k). The basic idea 
behind the calculation of Dij is to compare i-th polarization 9i with the 
n„ nearest neighbours 0^ of the j-th polarization 9j. Hence one possible 
coordinate invariant definition of Dij is 



k=l 

" COS [29i -2{9k + Afc^,)] (18) 

k=l 
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where Ok, k = l,...,n^ are the Uy nearest neighbours of the polarization 
9j, excluding 6j itself. Hence aU the Uy polarizations are compared to the 
polarization 9i after parallel transporting them to the i-th position. It is clear 
that if we replace the generalized product with the ordinary dot product in 
this equation, i.e. if we set A^^j = 0, then Dfj reduces to the matrix Dy 
that is used in the Andrews & Wasserman test with the normalization factor 
J\f in Eq. ^jset equal to unity. Dfj is the measure of alignment between the 
i-th vector and n„ nearest neighbours of the j-th vector. The statistic is then 
calculated in a manner similar to that used in the Andrews &: Wasserman 
test. We again arrange Dfp j = 1, ..,n, in ascending order and obtain the 
rank TZi of Dfj^^. The statistic is calculated using equation [T71 and is called 
ZP. 

Another somewhat more complicated coordinate invariant definition of 
Dij can be obtained by parallel transporting all the polarizations to the 
position of the polarization Vj and then comparing them at this point. The 
resulting expression is given by, 

= — E cos {[20 j -2(0^ + Ak^j)] - [20j -2{9i + A,^,)]} (19) 
k=i 

This can also be explicitly expressed in terms of the basic coordinate in- 
variant quantity, (Oi, Oj), but the precise expression is not needed here. The 
results obtained by using Dfj are similar to those obtained by using Dfj. 

3.1 Coordinate Invariance 

We next test explicitly whether the statistics defined in the previous section 
are indeed invariant under coordinate transformations. Consider a polariza- 
tion which is oriented at an angle and located at the position (a, 6) in a 
particular coordinate system. Next we consider a transformed coordinate 
frame such that its pole lies at the position (ap,Sp) in the original coordi- 
nates. Let the polarization be oriented at an angle On in the new system. 
The angle On can be expressed as (Hutsemekers 1998) 

^ N cos 5j) sm(aj) — a) 

tan {0 - On) = . a A \ V " ^^^^ 

sm Op cos — sm o cos dp cos(ap — a) 

The angular position of the polarization can also be evaluated in the new 
coordinate system by making an orthogonal transformation. The coordi- 
nates are transformed by first making a rotation about x-axis by an angle 
(— '5p+vr/2) and then a rotation about the new z-axis by an angle (— ap+7r/2). 
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The resulting statistics and for several choices of the pole position 
are given in Table 1, choosing the number of nearest neighbours = 20. 
The coordinate independence of the statistic is clear from this table. 



OLp 


8p 




ZP 





90 


0.267563 


2.069346 


100 


20 


0.267581 


2.069346 


25 


-80 


0.267569 


2.069346 


150 


55 


0.267586 


2.069346 


30 


-45 


0.267568 


2.069346 



Table 1: The statistics & for different choices of arbitrary poles 
{ap,6p) with riy = 20 nearest neighbours. 



4 Results 

The significance levels (S. L.) for the and statistic are given in Figures 
ISlini We plotted the logarithmic significance level of both the statistical tests 

and as a function of the number of nearest neighbours. The S. L. 
was computed by numerical simulations which compares the data statistic 
with that of a large number Ns of random samples. For most cases the choice 
Ns = 1000 was found to be sufficient. In several cases it was found necessary 
to increase Ns to 10000 or larger in order to get a reasonable estimate of 
S.L.. In computing Z^ we used the matrix Dfj given in Eqn. (|18|) . The 
matrix Df^, Eqn. gave similar results. The results are obtained both 

by using the radial distance r{z) given in Eqn. ()13|) as well as by setting 
r(z) = 1, which tests for a redshift independent alignment. We also try out 
several cuts based on the redshift z and the degree of polarization p. In 
Figs. I3|4I we show the significance level using the S'^, statistic for redshift 
dependent and redshift independent alignment respectively. In Figs. I5l6l we 
show the signifigance level using the Z^ statistic for redshift dependent and 
redshift independent alignment respectively. 

For an entire data set of 213 points we find that the log(S'.L.) < —2.5 
over a wide range of values using the test for redshift dependent 
alignment. The minimum log(S'.-L.) = —2.92 for = 32. Hence we find 
that there is evidence for large scale alignment as observed by Hutsemekers 
(1998). The ZP statistic (Fig. Ej) also shows alignment with log{S.L.) < -2 
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for > 32. The redshift independent test (Fig. |l]and|ni) does not show any 
evidence of ahgnment for the complete sample. We also examined several 
cuts on the data sample in order to determine whether the alignment arises 
from large or small degrees of polarization and redshifts. We find that low 
polarizations p < 2% show a very significant redshift dependent alignment 
as seen in Figures |21 and El The number of data points with this cut is 146. 
The log(5.L.) = —4.0 for = 24 and 32 using the statistic and reaches 
—5.0 for n„ = 28 with the statistic. The data with large polarization, 
p > 1%, however, shows no evidence of alignment with any choice of the 
statistic. The number of data points in this case is 147. The fact that low 
polarizations are significantly aligned suggests that the effect may arise due 
to propagation. 

We next compute the mean distance at which the alignment effect is 
most pronounced by directly testing for alignment over a certain distance, 
rather than using the number of nearest neighbours. In this case we test for 
alignment of a particular object with all the objects within a fixed distance 
of this object. Using the statistic with p < 2% we find that the S.L. takes 
its minimum value at r(z) ~ 0.22{2c/Hq), which corresponds to a distance 
scale of the order of a Gpc. 

We next examine the cuts on redshifts. In Fig. 13 we find that the 
sources at large redshifts z > 1 show a very strong alignment with the 
Sjy statistic for number of nearest neighbours larger than about 35. The 
number of objects in the data sample z > 1 are 115. For comparison we also 
consider a low redshift data sample, where the cut z < 1.3 is chosen such 
that it also contains approximately the same number of objects as the z > 1 
sample. The z < 1.3 infact also contains 115 objects. From Figs. |21andE]we 
find that the large redshift sample shows a strong alignment independent 
of whether we take the radial distance of the object into account or set 
it equal to unity. Hence the alignment of objects at large redshifts is not 
necessarily redshift dependent. The low redshift objects, however, do not 
show a strong alignment with this statistic. The results for the statistic 
(Figs. IS] and 121), however, show no alignment for large redshift sources and 
a weak alignment for low redshifts. The fact that the S'jj statistic shows 
a strong signal for alignment of large redshift objects and the Z^ shows 
no signal is easily understood. We find that the alignment of the large 
redshift sources occurs predominantly for relatively large number of nearest 
neighbours, and hence for large distances. The Z^ statistic dominantly 
tests only for local alignment. The can test both local and large scale 
alignment. The large redshift points seem to be aligned over very large 
distances. This is shown more clearly in Fig. which shows the log(S'.L.) 
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for very large values of Uy. The results for p < 2% cut are also shown 
for comparison. It is clear from this figure that the large redshift points 
show alignment over the entire sample. Hence we find that besides the 
redshift dependent alignment, which happens primarily for the objects with 
low polarization, the entire set of large redshift objects are aligned with one 
another. This effect is different from the redshift dependent alignment and 
was not noticed by Hutsemekers (1998). We emphasize that the redshift 
dependent alignment is seen dominantly for low polarizations as shown in 
fig. 01 The large redshift objects, however, show an alignment over the 
entire sky. 

In order to understand this very large scale alignment at large redshifts, 
z > 1, we make a scatter plot of the objects which show a significant align- 
ment along with those which do not. In this study we ignore the redshift 
dependence of the objects since the z > 1 set shows significant effect ir- 
respective of whether we include or ignore the redshift dependence. From 
fig. [7|we find that the S.L. is minimum at n„ = 38 and hence we choose 
this value of n„ for our study. The resulting scatter plot is shown in fig. |H1 
where the pluses and dots refer to the objects for which the dispersion mea- 
sure di < 0.25 and di > 0.25 respectively. We point out that the statistic 
Sjj = 0.307 in this case. We choose di = 0.25 as the cutoff in this figure since 
it shows the boundary between the aligned and non-aligned sources clearly. 
A choice of di = 0.3 as the cutoff also leads to a similar scatter graph, but 
with a larger overlap between the aligned and non-aligned regions. Further- 
more the figure does not change too much if we choose a larger value of 
n„. It is clear from figure |S1 that a large contribution to statistic is obtained 
from what is called the Al region by Hutsemekers (1998). This region is 
centered roughly at the Virgo supercluster and is delimited in redshift such 
that 1 < z < 2.3. However as can be seen from Fig. |H1 other regions also 
show significant alignment. 

In order to understand the nature of this very large scale alignment 
we first observe that the objects lying within the coordinate interval 10 < 
RA < 16 and —30 < Dec < 30 are aligned with one another. Similarly 
the objects in the region RA > 22, RA < 2 and —40 < Dec < 30 are also 
aligned over this entire patch. Hence we find that the data splits into two 
large patches such that in each patch most of the objects show alignment 
with one another. The mean polarization angles, i.e. the mean over the 
n„ = 38 nearest neighbours, in the interval 10 < RA < 16 are centered at 
the value of approximately 2.95 radians. In the region RA > 22, RA < 2, 
the means are centered around the value 2.35 radians. A parallel transport 
from RA = 0, Dec = —5.0 to RA = 13, Dec = leads to a shift in the 
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polarization angle by 0.64 radians, which when added to 2.35 leads to an 
angle very close to 2.95. Hence we find that even these two widely separated 
regions are correlated. This explains the significant alignment observed in 
the entire data sample corresponding to 2; > 1.0 as shown in fig. [7| 

5 Physical Explanation 

The observation of such a large scale alignment is quite surprising and is not 
easily explained in terms of conventional astrophysics. QSOs at such large 
distances from one another are unlikely to be correlated with one another 
and hence the effect most likely arises due to propagation. The fact that 
the alignment is most dominant in the data sample with small polarizations 
further supports this hypothesis. Polarizations can show alignment with one 
another if they arise dominantly due to extinction. For example the galactic 
magnetic fields result in a large scale alignment of the dust particles, which 
preferentially attenuate a particular polarization component of the electro- 
magnetic waves. Hence if the intrinsic source polarization is negligible, this 
will give rise to polarizations which are aligned over large distances. Galac- 
tic extinction, however, is unable to explain the redshift dependence of the 
effect. Furthermore it cannot explain why the large redshift points are cor- 
related over such large distances, whereas the low redshift data shows no 
such correlation. We next examine the possibility that the alignment arises 
due to supercluster extinction. We assume that the magnitude of the super- 
cluster extinction is large enough to cause such an alignment. The redshift 
dependent alignment seen in the low polarization sample {p < 2%) can then 
be explained if we assume the presence of a few very large superclusters of 
distance scale of order 1 Gpc. The cosmological scale alignment seen in the 
large redshift sample is, however, not easily explained by this mechanism 
since we do not expect the presence of a supercluster which covers almost 
the entire celestial sphere at redshift of 1. Hence the observations are not 
completely explained in terms of galactic or supergalactic extinction. 

Another physical phenomenon that can potentially explain the observa- 
tions is the pseudoscalar-photon mixing (Hutsemekers & Lamy 2001; Jain, 
Panda & Sarala 2002). A light pseudoscalar particle is predicted by many 
extensions of the standard model of particle physics (for e.g. see Peccei 

6 Quinn 1977; Mann & Moffat 1981; Sachs 1982; Will 1989; Ahluwaha & 
Goldman 1993; Ralston 1995). Its mixing with photons and its astrophysical 
consequences have also been studied by several authors (for e.g. see Sikivie 
1983; Maiani, Petronzio & Zavattini 1986; Harari & Sikivie 1992; Das, Jain 
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&; Mukherjee 2001). Its coupling to photons is bounded by observations of 
SN87A to be c/ < IQ-^^ GcV~^ (Brockway, Carlson & Raffclt 1996; Grifols, 
Masso & Toldra 1996; Raffelt 1999; Rosenberg & van Bibber 2000). Such a 
particle has also be invoked to explain the observed dimming of supernovae 
at large redshifts even if the expansion rate of the universe is not acceler- 
ating (Csaki, Kaloper k, Terning 2002). The pseudoscalar particle decays 
into photons in the presence of background magnetic field 5, such that the 
photon produced is polarized parallel to the transverse component of B, 
which we denote as Bt- Hence the electromagnetic wave produced by this 
mechanism is polarized parallel to Bt- Similarly a photon polarized parallel 
to B can decay into the pseudoscalar during propagation, which leads to a 
wave polarized perpendicular to Bt- 

The existence of a hypothetical pseudoscalar can explain the observed 
alignment as follows. Wc first assume that the redshift dependent effect seen 
in low polarization objects is explained by the presence of a few superclusters 
of length scales of order 1 Gpc. As the electromagnetic waves pass through 
these superclusters the polarizations get aligned either due to extinction or 
due to pscudoscalar-photon mixing. Wc further assume that the large red- 
shift objects, which show alignment over the entire sky, emit a large flux of 
the pseudoscalar particles whereas the flux emitted by low redshift objects 
is relatively small. This is reasonable since the large redshift objects are 
in general very active and have high temperatures. As these QSOs evolve 
they are less active and the pseudoscalar flux becomes negligible, which ex- 
plains why the low redshift objects behave differently. As the pseudoscalars 
propagate through our local supercluster they decay into photons which are 
polarized parallel to the supercluster magnetic field. As discussed by Jain 
et al. (2002) the decay probability of such a pseudoscalar is of order unity 
with the current limits on the pscudoscalar-photon coupling. In obtaining 
this estimate we took the Virgo supercluster parameters for the magnetic 
field [B ^ 1/LtG) (Vallee 1990) and plasma density {ue ~ 10"^ cm'^^). We 
propose that a weaker magnetic field, B ~ 0.1/iG, might be associated with 
the entire supercluster, and may be responsible for the observed alignment 
of large redshift objects over very large angular separations. We point out 
that the decay probability of pseudoscalars into photons 7, P^^^ B^/ri^ 
(Carlson & Garretson 1994; Jain et al. 2002) and hence even if the magnetic 
field is an order of magnitude smaller than that observed in the Virgo su- 
percluster, P(f)^^ can be large as long as Ue also decreases proportionately. 
Hence this phenomenon explains the alignment seen in this sample over very 
large angles. 

It is clearly very important to further test the proposal that the align- 
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merit effect is due to pseudoscalar-pfioton mixing. Jain et al. (2002) liave 
computed all the Stokes parameters of the electromagnetic wave in the pres- 
ence of such a pseudoscalar particle. They find that the spectral dependence 
of the linear polarization, the circular polarization and the orientation angle 
of the linear polarization are all closely correlated with one another due to 
this mixing. This correlation can be studied by making further observations, 
which can establish or rule out this explanation. 

In our discussion we have assumed an isotropic and homogeneous uni- 
verse. We have not considered the possibility of a cosmological scale mag- 
netic field or a Lorentz violating interaction (Nodland & Ralston 1997; Car- 
roll, Field & Jackiw 1990), which might also explain these observations. 
A large scale dipole anisotropy is also observed in the radio polarizations 
from distant AGNs (Jain & Ralston 1999). This anisotropy was found to 
be independent of radial distances and might be explained in terms of some 
local effect. The axis of this dipole anisotropy is found to be pointing ap- 
proximately opposite to the center of our local supercluster. Since this is a 
dipole anisotropy the effect on polarizations is large both in the direction 
and opposite to the axis. Hence the radio anisotropy displays an intriguing 
relationship to the large scale alignment seen in optical polarizations in the 
direction of the supercluster center. This might be a hint for a common 
origin of these two effects. The fundamental origin of the radio anisotropy 
is not known. Jain et al. (2002) pointed out that the existence of a light 
pseudoscalar can also explain this effect if the AGNs emit a large flux of 
pseudoscalars at radio frequencies. Although this explanation is disfavored 
due to the very large pseudoscalar flux required from the distant AGNs, it 
does explain the correlation between the radio dipole axis and the observed 
large scale alignment in optical polarizations. 

6 Conclusions 

In conclusion we have developed coordinate invariant statistical procedures 
in order to test for the large scale alignment of optical polarizations. We 
applied these tests to a data sample of QSOs compiled by Hutsemekers 
(1998) and Hutsemekers & Lamy (2001) and find that the polarizations 
show statistically significant alignment over very large distances. We find 
that the alignment is redshift dependent and is seen dominantly for the 
data sample with low polarizations {p < 2%). We also find that the large 
redshift, z > 1, sample shows a very large scale alignment. Infact the 
polarizations contained in almost the entire data sample at z > 1 seem 
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to be correlated with one another. We find that galactic or supercluster 
extinction is unlikely to provide an explanation for these observations. We 
also argue that the existence of a hypothetical pseudoscalar particle might 
provide an explanation for the alignment effect. 
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Figure 3: The logarithmic significance level, log(5.L.), as a function of 
the number of nearest neighbours n„ using the statistic S^. The nearest 
neighbours are obtained by taking into account the radial distance of the 
source and hence this tests for redshift dependent alignment. The black 
curve corresponds to the entire data set. The short dashed, dash-dotted, 
long dashed and dotted curves correspond to the cuts p < 2%, p > 1%, 
z > 1.0 and z < 1.3 respectively. 
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Figure 4: The logarithmic significance level, log(5.L.), as a function of 
the number of nearest neighbours n„ using the statistic S^. The nearest 
neighbours are obtained without taking into account the radial distance of 
the source and hence this tests for redshift independent alignment. The 
black curve corresponds to the entire data set. The short dashed, dash- 
dotted, long dashed and dotted curves correspond to the cuts p < 2%, 
p > 1%, z > 1.0 and z < 1.3 respectively. 
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Figure 5: The logarithmic significance level, log(S'.L.), as a function of 
the number of nearest neighbours using the statistic Z^. The nearest 
neighbours arc obtained by taking into account the radial distance of the 
source and hence this tests for redshift dependent alignment. The black 
curve corresponds to the entire data set. The short dashed, dash-dotted, 
long dashed and dotted curves correspond to the cuts p < 2%, p > 1%, 
z > 1.0 and z < 1.3 respectively. 
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Figure 6: The logarithmic significance level, log(S'.L.), as a function of 
the number of nearest neighbours using the statistic Z^. The nearest 
neighbours are obtained without taking into account the radial distance of 
the source and hence this tests for redshift independent alignment. The 
black curve corresponds to the entire data set. The short dashed, dash- 
dotted, long dashed and dotted curves correspond to the cuts p < 2%, 
p > 1%, z > 1.0 and z < 1.3 respectively. 
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Figure 7: The logarithmic significance level, log(S'.L.), as a function of 
the number of nearest neighbours using the statistic for the cuts 
z > 1 (solid curve) and p < 2% (dashed curve). The nearest neighbours are 
obtained without taking into account the radial distance of the source and 
hence this tests for redshift independent alignment. Results are shown for 
very large values of and this tests for alignment over very large distances. 
The total number of points in the sets z > 1 and p < 2% are 115 and 146 
respectively. 
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Figure 8: Scatter plot of the objects which show significant ahgnment di > 
0.25 (dots) and those which do not show ahgnment di < 0.25 (pluses). Here 
di is a measure of the dispersion as defined in Eq. 1141 The figure shows 
the data sample with the cut z > 1 with the number of nearest neighbours 
= 38. 
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